Exclusive processes in electron - ion collisions 
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The exclusive processes in electron-ion (eA) interactions are an important tool to investigate 
the QCD dynamics at high energies as they are in general driven by the gluon content of the 
target which is strongly subject to parton saturation effects. In this paper we compute the cross 
sections for the exclusive vector meson production as well as the deeply virtual Compton scattering 
(DVCS) relying on the color dipole approach and considering the numerical solution of the Balitsky- 
Kovchegov equation including running coupling corrections (rcBK). The production cross sections 
obtained with the rcBK solution and bCGC parametrization are very similar, the former being 
slightly larger. 

PACS numbers: 12.38.-t, 24.85. +p, 25.30.-c 

Exclusive processes in deep inelastic scattering (DIS) have appeared as key reactions to trigger the generic mecha- 
nism of diffractive scattering (For a recent review see, e.g. [l|). In particular, diffractive vector meson production and 
deeply virtual Compton scattering (DVCS) have been extensively studied at HERA and provide a valuable probe of 
the QCD dynamics at high energies. These processes are driven by the gluon content of the target (proton or nucleus) 
which is strongly subject to parton saturation effects as well as to nuclear shadowing corrections when one considers 
scattering on nuclei (See e.g. Q). In particular, the cross sections for exclusive processes in DIS are proportional to 
the square of the scattering amplitude, which makes them strongly sensitive to the underlying QCD dynamics. In 
a recent paper Q we have estimated the coherent and incoherent cross sections for exclusive p and J/\& production 
considering the color dipole approach and phenomenological saturation models which describe the scarce F A data 
as well as the HERA data. Our results demonstrated that the coherent production of vector mesons is dominant, 
with a small contribution coming from incoherent processes. Moreover, our results indicate that the experimental 
study of these processes is feasible in future electron - ion collider, as e.g. the eRHIC [lj or LHeC 5]. In this paper 
we complement our previous analysis including <f> production and extending our study to the nuclear DVCS (See 
also @). Moreover, we review our results for p and Jf^l production making use of the numerical solution of the 
Balitsky-Kovchegov equation 0, Q including running coupling corrections [^4ll| in order to estimate the contribution 
of the saturation physics to exclusive processes. Our main motivation is associated to the fact that the improved BK 
equation has been shown to be really successful when applied to the description of the ep HERA data on inclusive 
and diffractive proton structure function [l2l - fl4j |. as well as on exclusive processes [l5| and on the forward hadron 
spectra in pp and dA collisions [ij], [H| ■ 

Let us start presenting a brief review of exclusive processes in electron - ion collisions (For details see [H, [l7| ) ■ In 
the color dipole approach the exclusive production j*A — > EY (E = p, <f), J/fy or 7) in electron-nucleus interactions 
at high energies can be factorized in terms of the fluctuation of the virtual photon into a qq color dipole, the dipole- 
nucleus scattering by a color singlet exchange and the recombination into the exclusive final state E. This process 
is characterized by a rapidity gap in the final state. If the nucleus scatters elastically, Y = A, the process is called 
coherent production and the corresponding integrated cross section is given in the high energy regime (large coherence 
length: Z c > R A ) by |,[l7| 



eoh (<fA->EA)= d 2 b(M A (x,r,b)) 2 



(1) 
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where 

(AO = / d2r I dzV* E (r 7 z)Ar A (x,r,b)t> r (r,z,Q 2 ) (2) 



and Af(x, r, b) is the forward dipole-target scattering amplitude for a dipole with size r and impact parameter b which 
encodes all the information about the hadronic scattering, and thus about the non-linear and quantum effects in the 
hadron wave function. On the other hand, if the nucleus scatters inelastically, i.e. breaks up (Y = X), the process 
is denoted incoherent production. In this case one sums over all final states of the target nucleus, except those that 
contain particle production. The t slope is the same as in the case of a nucleon target. Therefore we have: 

a- ( 7 *A -> EX) = |Jm - 4(M = ° )|2 (3) 

1d7t Be 

where at high energies (l c >• R A ) [13] : 

\lmA\ 2 = J d 2 bT A {b) (a dp exp[~a dp T A (p)]\ (4) 

and a dp is the dipole-proton cross section, which in the eikonal approximation it is given by: 

a dp (x,r) = 2 J d 2 bAfP{x,r,b) . (5) 

In the incoherent case, the qq pair attenuates with a constant absorption cross section, as in the Glauber model, 
except that the whole exponential is averaged rather than just the cross section in the exponent. As discussed in 
the coherent and incoherent cross sections depend differently on t. At small-t (~tR A /3 <C 1) coherent production 
dominates, with the signature being a sharp forward diffraction peak. On the other hand, incoherent production will 
dominate at large-i (—tR A /3 ^> 1), with the i-dependence being to a good accuracy the same as in the production 
off free nucleons. 

In the Eqs. ([2]) and (HJ the functions ^(z, r) and & (z, r) are the light-cone wavefunctions of the photon and 
the exclusive final state, respectively. The variable r defines the relative transverse separation of the pair (dipole) 
and z (1 — z) is the longitudinal momentum fraction of the quark (antiquark). In the dipole formalism, the light-cone 
wavefunctions ^{z, r) in the mixed representation (r, z) are obtained through a two dimensional Fourier transform of 
the momentum space light-cone wavefunctions ^(z, k). The photon wavefunctions are well known in literature [l8j]. 
For the meson wavefunction, we have considered the Gauss-LC model of Ref. [18] . The motivation for this choice is its 
simplicity and the fact that the results are not very sensitive to differences between the models analyzed in [18j ]. We 
choose the quark masses to be m u ^d, s = 0.14 GeV and m c = 1.4 GeV. The parameters for the meson wavefunction can 
be found in Ref. [18|. In the DVCS case, as one has a real photon at the final state, only the transversely polarized 
overlap function contributes to the cross section. Summed over the quark helicities, for a given quark flavour / it is 
given by [H|, 

N a e 2 

- C 27T "° f { [z 2 + z 2 } e x K x {e X T)e 2 K x [e^) + m}K ( £l r)K (s 2 r)} , (6) 

where we have defined the quantities e\ 2 = zzQ\ 2 + m 2 and z = (1 — z). Accordingly, the photon virtualities are 
Q\ = Q 2 (incoming virtual photon) and Q 2 . = (outgoing real photon). 

In order to estimate the coherent production in eA collisions we need to specify the forward dipole - nucleus 
scattering amplitude, Af A (x, r, b). Following Q we will use in our calculations the model proposed in Ref. [2(J, which 
describes the current experimental data on the nuclear structure function as well as includes the impact parameter 
dependence in the dipole nucleus cross section. In this model the forward dipole-nucleus amplitude is given by 



JV A (x,r,b) = 1 -exp 



~a dp (x,r 2 )T A (b) 



(7) 



where o dv is the dipole-proton cross section and T A (b) is the nuclear profile function, which is obtained from a 3- 
parameter Fermi distribution for the nuclear density normalized to A. The above equation sums up all the multiple 
elastic rescattering diagrams of the qq pair and is justified for large coherence length, where the transverse separation 
r of partons in the multiparton Fock state of the photon becomes a conserved quantity, i.e. the size of the pair r 
becomes eigenvalue of the scattering matrix. 

In our approach the coherent [Eq. (JTJ] and incoherent [Eq. cross sections can be calculated in terms of the 
dipole-proton cross section or the forward dipole-proton scattering amplitude [See Eq. ([5])]. which is a solution of the 
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Q 2 = 1 GeV 2 
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FIG. 1: (Color online) Energy dependence of the coherent cross section at different final states and Q 2 — 1 GeV 2 



BK equation. As the leading order solution of the BK equation was not able to describe the HERA data, in Ref. Q 
we have used the GBW [2l[ and bCGC [l8| parametrizations for Af p as input in our calculations. However, in the last 
years the next-to-leading order corrections to the BK equation were calculated [H-EH • Such calculation allows one to 
estimate the soft gluon emission and running coupling corrections to the evolution kernel. The authors have verified 
that the dominant contributions come from the running coupling corrections, which allow us to determine the scale of 
the running coupling in the kernel. The solution of the improved BK equation was studied in detail in Refs. (ToT l22l|. 
Basically, one finds that the running of the coupling reduces the speed of the evolution to values compatible with 
experimental data. In (l2j a global analysis of the small x data for the proton structure function using the improved 
BK equation was performed (See also Ref. (H)]). In contrast to the BK equation at leading logarithmic a s ln(l/x) 
approximation, which fails to describe the HERA data, the inclusion of running coupling effects in the evolution 
renders the BK equation compatible with them (See also [141-1 1610 . It is important to emphasize that the impact 
parameter dependence was not taken into account in Ref. [12J, the normalization of the dipole cross section was 
fitted to data and two distinct initial conditions, inspired by the Golec Biernat-Wusthoff (GBW) [2l[ and McLerran- 
Venugopalan (MV) [23[ models, were considered. The predictions resulted to be almost independent of the initial 
conditions and, besides, it was observed that it is impossible to describe the experimental data using only the linear 
limit of the BK equation. The parametrizations obtained in (l2j were very successful in reproducing DIS data but it 
remains to seen whether they can also be used to describe data from RHIC. Other parametrizations of dipole cross 
sections had to be slightly modified in order to account for RHIC data [24| . 

In what follows we calculate the exclusive observables using as input in our calculations the solution of the running 
coupling Balitsky-Kovchegov (rcBK) evolution equation. We make use of the public- use code available in [25[. In 
numerical calculations we have considered the GBW initial condition for the evolution. Furthermore, we compare 
the rcBK predictions with those from the bCGC model used in our previous calculation Q. Moreover, in order to 
calculate the incoherent cross section for vector meson production we will use the following parametrization of the 
diffractive slope 



B v (Q 2 ) = 0.60 



14 



(Q 2 



M 2 ) - 26 



1 



(8) 



obtained from a fit to experimental data referred in Ref. [26J . In the DVCS case we take the experimental parametriza- 
tion H3, B (Q 2 ) = a[l - blog(Q 2 /Q 2 )}, with a = 6.98 ± 0.54 GeV 2 , b = 0.12 ± 0.03 and Q 2 = 2 GeV 2 . 

In Fig. [T]we show the coherent production cross section as a function of the photon-target c.m.s energy, W, for a 
fixed photon virtuality Q 2 — 1 GeV 2 . Each one of the panels shows the results obtained for one specific final state. 
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(Color online) Dependence on the photon virtuality of the coherent cross section for different final states and W = 500 



In each single figure the two upper (lower) curves show the results for a Pb (Ca) target. In all figures the dashed 
(solid) lines are obtained with the bCGC (rcBK) dipole cross section. Fig. [2] shows the same cross sections, this time 
as a function of Q 2 for a fixed energy, W — 500 GeV. Figs. [3] and 0] are the exact analogues (of Figs. [T]and[2]) for the 
corresponding incoherent cross sections. 

The curves in the figures have the merit of being the first concrete predictions made for these processes with the 
help of the recently obtained rcBK dipole cross section. They present some features which are expected and some 
other features which could not have been anticipated without a quantitative calculation. In first place we observe, as it 
should be, that all cross sections grow with W and fall with Q 2 . The first feature is related solely to the nature of the 
dipole cross section, which grows with the energy, whereas the second feature comes from the dipole wave functions. 
We can also see from the figures that, at least for the two cases considered (bCGC and rcBK), the production cross 
sections are not very strongly dependent on the choice of the dipole cross section. 

At low Q 2 and low W the bCGC and rcBK production cross sections are indistinguishable one from the other 
because the dipole cross sections tend to coincide. These latter have been tuned to fit DIS data, which are taken in 
this kincmatical region. Another expected feature is the observed decrease of the cross sections with increasing vector 
meson masses, which comes from the wave functions. 

Differences are expected to appear at higher energies, where we enter the lower x (extrapolation) region. In all 
cases we see that the results obtained with the rcBK cross section are larger than those obtained with the bCGC one. 
This is related to the fact that the numerical solutions of the BK equation tend to reach later the unitarity limit [28| . 
In the first estimates, with a fixed coupling, the solutions of the BK equation would saturate too fast. In subsequent 
studies it was found that running coupling corrections to the BK kernel could bring the evolution speed down to 
values compatible with those extracted from data, but still larger than those found in other parametrizations, such 
as the bCGC one. Due to this fact, the results obtained with the rcBK dipole cross section grow faster with energy 
than those obtained with the bCGC one. 

A curious feature in the figures is that the differences between bCGC and rcBK are larger for heavier vector mesons. 
This can be understood looking carefully at the integrand of (JTJ, which is the product of the wave functions, containing 
information about the masses, and the dipole cross section. As a function of the dipole size r the difference between 
bCGC and rcBK is mostly in the low to intermediate r region, where the bCGC is always below the rcBK dipole 
cross section. At large r the two cross sections are close to each other. The overlap function, i.e., given by 

© (with the inclusion of the longitudinally polarized overlap function), has peaks at different locations. The p is a 
larger state and its overlap function peaks at much larger values than the J ftp overlap function. In this way it gives 
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FIG. 3: (Color online) Energy dependence of the incoherent cross section for different final states and Q 2 = 1 GeV 2 . 

a stronger weight to larger r where the differences between bCGC and rcBK are smaller. The same thing happens 
to the <j). On the other hand, the J/ip overlap function peaks at smaller r where the dipole cross sections are more 
different from each other. A similar behavior is verified in the DVCS case. 

As shown in [28], at increasing values of the energy W (and thus of smaller of x) the difference between rcBK 
and bCGC moves to smaller values of r, a region which is suppressed by the overlap functions of the p and <j>. This 
explains why the p (and also the <j>) production cross sections are the almost the same for rcBK and bCGC dipole 
cross sections for all the energies considered, as it can be seen on the right side of Figs. 1 and 3. 

In Figs. 2 and 4 we would expect to see a convergence of the curves for higher values of Q 2 . In this region the 
dipoles are small and all cross sections should approach the color transparency regime. In fact, a difference between 
them persists even at large Q 2 because in the expressions used here there is no DGLAP evolution, which would bring 
the dipole cross sections together. 

As a summary, we presented a systematic analysis of exclusive production in small- x deep inelastic electron - ion 
scattering in terms of the non-linear QCD dynamics. This was the first calculation (of these observables) using the 
solution of the BK equation improved with running coupling corrections. In this work we obtain predictions for the 
exclusive production of vector mesons and DVCS. Our analysis confirms the dominance of the coherent production 
with a small contribution coming from incoherent processes, a result previously found in Q. Our main result is 
that the BK evolution equation implies larger cross sections for exclusive processes than the phenomenological model 
proposed in [l8j], the so called bCGC model. Our predictions for both vector meson and DVCS production are relevant 
for the physics programs of the future experiments eRHIC and LHeC. 
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